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Abstract 

The paper is devoted to retrieval of the first order phase transition signal in the 
inelastic collisions. The primary intent is to show that the experimentally observ- 
(yQ • able signal exist iff the multiplicity is sufficiently large. We discuss corresponding 

^D , phenomenology from the point of view of experiment. 
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The question of possibility to observe first order phase transition in the hadron and 
heavy ion collisions is discussed widely at present time [T]. The aim of paper is discuss 
the phenomenology of that problem. 

The first order phase transition in statistics appears in the result of creation of local 
^ ■ critical fluctuations (e.g. bubbles of vapor if the (liquid— >gas) transition is considered), 

in contrast to the second order phase transition where the whole system undergo the 
0-(, transition. The dimension of such fluctuations increase, if they are under-critical, and 

r-| I the whole system in result undergo the transition. 

It is hard to imagine that exactly such picture appears in the hadron or ion inelastic 
CN ■ collisions. The point is that in the event-by-event measurement the dimension of bub- 

-^ . bles, if they exist, may be both smaller or larger critical one. Therefore, in the best case, 

OA ! we observe the mixture of two-phase medium and the question how one may increase 

on experiment the weight of events with under-critical bubbles is the flrst problem. It is 
evident that the density, usually used in statistics, can not be introduced as the "order 
Q ■ parameter" since only the result of particles production process is observed. 

00 , In the paper [2] the "chemical potential", fi{n, s), was offered as the "order parame- 

ter" . It is the work which is necessary for creation of one particle and it was shown that 
if the role of under-critical babbles dominate then fi{n, s) must decrease with number 
/\ I of produced particles n. In another words, it was offered to deflne the boiling, i.e. cre- 

c^ ■ ation of under-critical bubbles, through intensity of the process of particles production 

(evaporation). 

This is the main idea why the very high multiplicity (VHM) processes were consid- 
ered. It is not important what additional criterium is used searching flrst order phase 
transition. In any case we must consider the VHM domain to have intensive production 
of particles. In addition, the kinetic degrees of freedom must be suppressed in the VHM 
region. 

Therefore, the special attention will be given to the VHM processes. Then exist 
approximation p]: 

/i(n, s) ~ '- — lno-„(s) (1) 

n 

in this multiplicity region. Here T is the mean energy, including mass, of produced par- 
ticles, i.e. T is associated with temperature, and (j„ is the normalized to unite multiple 
production cross section which can be considered in the VHM region as the "partition 
function" of the equilibrium system. The equilibrium condition will be deflned later, 
see inequality ( TT7I) . deflnition (IT8l) and [3] where the detailed explanation was given. See 
also the footnote 5. 
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Continuing the analogy with thermodynamics one can say that (— Tlno"„)/n is the 
Gibbs free energy per one particle. Then /x can be interpreted as the "chemical potential" 
measured with help of observed free particleqj. 

The definition ([T]) is quiet general. It can be used both for hadron-hadron and ion-ion 
collisions, both for low and high energies, both for "boiling" media of coloured partons 
and colorless hadrons. Definition ([1]) is model free and operates only with "external" 
directly measurable parameters. The single indispensable condition: we work in the 
VHM region. 

We will discuss in the paper the chance of experimental measurement of fi{n, s) 
defined by ([T]), what kind uncertainties hides it from experimental point of view noting 
the the cross sections in VHM domain are small. The correction to ([1]) are not essential 
in the VHM region but nevertheless the field-theoretical definition of chemical potential 
in using Wigner functions formalism [4j| will be published, see also [3]. 

II. 

It is necessary to remind main steps toward ([T]) to understand hidden phenomeno- 
logical uncertainties. The starting point [2J was the generating function 

oo 

p{z, s) = ^ 2;"cr„(s), p(l, s) = 1, an = at n> n^ax = Vs/m, (2) 

n=0 

where m is the hadron mass. One may use inverse Mellin transformation: 

to find o"„ if p{z, s) is known. One may calculate integral ([3]) by saddle point method. 
The equation (of state): 

d 

n = z—\iip{z,s) (4) 

defines mostly essential value z = z{n, s). Therefore, only 

^ ^ ^max zyn^jinx, S) yo) 

have the physical meaning. 

One may write p{z, s) in the form: 

p(z,s)=exp|^A(5)[, (6) 

where the "Mayer group coefficient" 6/ can be expressed through correlators Ck{s): 

Let us assume now that we have Poisson distribution, i.e. if in the sum: 

lnp(2;, s) = Y^ ^-^j Ck{s) (7) 

k 

one may leave first term, then it is easy to see that 

z{n,s) =n/ci{s), ci{s) = n{s), (8) 



^Notice that one may consider n as the muhipUcity in the experimentaUy observable range of phase 
space. 



are essential and in the VHM region: 

Tl 

lna„(s) = -nln^— (1 + 0(1/ Inn)) = -nlnz(n, s)(l + 0(l/lnn)). (9) 

Therefore, in considered case with c^ = 0, A; > 1, exist following asymptotic estimation 
for n » 1: 

IncTrt — — rilnz(n, s), (10) 

i.e. an is defined in VHM region mainly by the solution of Eq.(j4]) and the correction 
can not change this conclusion. It will be shown that that kind estimation is hold for 
arbitrary asymptotics of cr„. The definition ([1]) based on this observation. 

If we understand an as the "partition function" in the VHM region then z is the 
activity usually introduced in statistical physics. Correspondingly the chemical potential 
H is defined trough z: 

fx = Tlnz. (11) 

Combining this definition with estimation ( ITOl) we define an through /x. But, if this 
estimation does not depend from asymptotics of o"„ over n, i.e. if it has general meaning, 
then it can be used for definition of /x(n, s) through o"„(s) and T{n, s) at n » 1. Just 
this formal idea is realized in ([1]): it can be shown in Sec. HI that ([1]) is correct at the 
asymptotical value of n. 

III. 

Now we will make the important step. To put in a good order our intuition it 
is useful to consider p{z, s) as the nontrivial function of z. In statistical physics the 
thermodynamical limit is considered for this purpose. In our case the finiteness of energy 
■\/s and of the hadron mass m put obstacles on this way since the system of produced 
particles necessarily belongs to the energy-momentum surfacqj. But we can continue 
theoretically an to the range n > rimax and consider p{z, s) as the nontrivial function of 
z. This step hides the assumption that nothing new appear at n > rimax, i-e. the VHM 
interval -n << n < rimax is sufficiently wide to represent main physical processes. 

Let us consider the analog generating function which has the first n < rimax coeffi- 
cient of expansion over z equal to an and higher coefficients for n > rimax are deduced 
from continuation of theoretical value of cr„ to n > rimax- Then the inverse Mellin 
transformation ([3]) gives a good estimation of cr„ through this generating function if the 
fluctuations near zin, s) are Gaussian or, it is the same, if 



\2n — z^d'^\np{z, s)/dz^ 



\n + z'^d'^ \np{z, s)/9z^|^/^ 



« 1. (12) 

-z{n,s) 



Notice that if the estimation (ITOl) is generally rightful then one can easily find that l.h.s. 
of flT2|) is ~ l/n^/"^. Therefore, one may consider p{z^ s) as the nontrivial function of z 
considering z{n, s) < Zmax \in«n < rimax- 

Then it is easily deduce that the asymptotics of crn(s) is defined by the leftmost sin- 
gularity, Zc-, of function p(^, s) since, as it follows from Eq.(jll), the singularity "attracts" 
the solution z[ri^s) in the VHM region. In result we may classify asymptotics of an in 
the VHM region if ([12]) is hold. 

Our problem is reduced to the definition of possible location of leftmost singularity of 
p{z, s) over 2; > Qj. It must be stressed that the character of singularity is not important 



■^It must be noted that the canonical therniodynaniic system belongs to the energy-momentum shell 
because of the energy exchange, i.e. interaction, with thermostat. The width of the shell is defined by 
the temperature. But in particle physics there is no thermostat and the physical system completely 
belongs to the energy momentum surface. 

"^The singularities in complex z plane will not be considered since they lead only to oscillations in 
multiplicity distribution. 



for definition of ^{n,s) in the VHM region at least with 0(1/ Inn) accuracy. One may 
consider only three possibility at n — i> oo: (I) z{n, s) ^ Zc = 1; (II) z{n, s) ^ z^, 1 < 
Zc < cxd; (III) z{n, s) ^> Zc = cxd. The structure of complex z plane is much more 
complicate but for our purpose the above described picture is sufficient. 

Correspondingly one may consider only three type of asymptotics in the VHM region: 
(I) o"„ > 0(e~"); (II) a„ = 0(e~"). Such asymptotics is typical for hard processes 
with large transverse momenta, like for jets [3J; (III) o"„ < 0(6""). That asymptotic 
behavior is typical for multiperipheral-like kinematics, where the longitudinal momenta 
of produced particles are noticeably higher than the transverse ones ^ . 

Therefore the case (I) is the best candidate for phase transition since in this case the 
cross sections are comparatively large in the VHM region, i.e. particles "intensively" 
produced in that case. Notice that if (I) is not realized in nature then the (II) kind 
processes would dominate in the VHM region. 

Let us consider now the estimation ([T]). It follows from ^ that, up to the preexpo- 
nential factor, 

lno"(ra, s) ~ —n\Yiz{n, s) + \n.p{z{n, s), s). (13) 

We want to show that, in a vide range of n from VHM region, 

nln2;(n, s) ~ lnp(2;(n, s), s). (14) 

Let as consider now the mostly characteristic examples. 

(I) Singularity at z = 1. The physical meaning of singularity at z = 1 may be 
illustrated by the droplet model [6]. The Mayer's group coefficient, see ([6]), for cluster 
from / particle is 6z(/3) ~ exp{— /?rZ*-'^~^'*''^}, where r/*-''"^^''^, / >> 1, is the surface tension 
energy, d is the dimension. Therefore, ii d > 1 the series over / in (Q diverges at z = 1. 

This case was considered in [2] in details. In the used lattice gas approximation 
hiz{n) ~ n~^ and lncr„ ^ — n~^ = — nln2;(n)(l + 0{l/n)). Notice that the simplest 
droplet model predicts unphysical asymptotics: (T„ —>■ const in the VHM region. 

(II) Singularity at 1 < Zc < oo. Let us consider one jet contribution: In p{z,s) = 
— 7ln(l — nj{s){z — 1)). In this case z{n,s) = Zc{l — 'j/n), n » 7, and lncr„ = 
—nlnz{n, s)(l + 0{\nn/n)). 

(HI) Singularity at z = 00. For k Pomeron exchange: \np{z,s) = Ck{s){z — 1)^. In 
this case z{n) = {n/kckY^'^ >> 1 and lncr„ ~ —n\nz{n)(l + 0(l/lnn)). 

One can conclude: 

(i) The definition ([T]) in the VHM region is rightful since the correction falls down 
with n. On this stage we can give only the qualitative estimation of corrections. Nev- 
ertheless (II]) gives the correct n dependence in the VHM region. 

(ii) Activity z{n, s) tends to Zc from the right in the case (I) and from the left if we 
have the case (II) or (HI). 

(iii) The accuracy of estimation of the chemical potential ([T]) increase from (HI) to 

(I)- 
IV. 

The temperature T is the next problem. The temperature is introduced usually using 
Kubo-Martin-Schwinger (KMS) periodic boundary conditions. But this way assumes 
from the very beginning that the system (a) is equilibrium and (b) is surrounded by 
thermostat through which the temperature is determined. The first condition (a) we 
take as the simplification which gives the equilibrium state. 

The second one (b) is the problem since there is no thermostat in particle physics. 
For this reason we introduce the temperature as the Lagrange multiplier f3 = 1/T 
of energy conservation law [5]. In such approach the condition that the system is in 



equilibrium with thermostat replaced by the condition that the fluctuations in vicinity 
of P are Gaussian. 

The interesting for us p{z, s) we define through inverse Laplace transform of p{z, (3): 

p(^,^)=/^^e^^^p(;.,/?). (15) 

It is known that if the interaction radii is finite, i.e. the hadron mass is finite, then the 
equation (of state): 

^s = -^p{z,(3) (16) 

have real positive solution P{n, s) at z = z{n, s). We will assume that the fluctuations 
near f3{n, s) are Gaussian. This means that the inequality [3]: 

\dHnp{z,(3)/d(3^\ 



\dHnp{z,(3)/d(3^\^/^ 



« 1 (17) 

z=z(n,s),l3=l3{n,s) 



is satisfled. Therefore, we prepare the formalism to find "thermodynamic" description of 
the processes of particle production assuming that this S'-matrix condition of equilibrium 
dUD is holci. 

I want to underline that our thermal equilibrium condition 017p have absolute mean- 
ing: if it is not satisfied then /3(n, s) loses every sense since the expansion in vicinity 
of /3(n, s) leads to the asymptotic series. In this case only the dynamical description of 
^'-matrix can be used. 

It is not hard to see [Sj that 

gl I 



d/3 



\nR{z,p)\^=;,(^ri,s),/3=l3{n,s) =< Y\.iVi- < V >) >n,s 



i=l 



is the /-point energy correlator, where < ... >n,s nieans averaging over all events with 
given multiplicity and energy. Therefore 0171) means "relaxation of /-point correlations", 
/ > 2, measured in units of the dispersion of energy fluctuations, 1 = 2. One can note 
here the difference of our definition of thermal equilibrium from thermodynamical one 

m 

V. 

We may conclude that: 

(i) The definition of chemical potential ([T]) was discussed. This important observable 
can be measured on the experiment directly. Chemical potential, p{n, s), must decrease 
in the VHM region if the first order phase transition occur, case (I), and it rise in 
opposite case, see (II) and (III), see Sec. III. 

(ii) We are forced to assume that the energy and the multiplicity are sufficiently 
large, i.e. the experimental value 2;'^^^(n, s) is sufficiently close to z^ = 1. In opposite 
case the leading leftmost singularity over z would not be "seen" on experiment and the 
production processes constitutes from the complicated mixture of subprocesses. 

(iii) The cross section an falls down rapidly with n and for this reason the VHM 
events are hardly observable. One may avoid this problem considering the finite en- 
ergy heavy ion collisions as the most candidates of processes described by methods of 
thermodynamics and Zc is easier "reachable" in this case. 



^Introduction of /3(n, s) allows to describe the system of large number of degrees of freedom in terms 
of single parameter /3(n, s), i.e. it is nothing but the useful trick. It is no way for this reason to identify 
entirely l//3(n, s) with thermodynamic temperature where it has the self-contained physical sense. 
Nevertheless path-integral representation of p{P, z) defined from S'-matrix coincides with Feynnian- 
Kac representation of grand partition function [3] if PT|) is hold. It must be noted also that the energy 
spectrum of produced particle in this case have Boltzmann form, e~^^ . 



(iv) One can define z{n,s) also directly from Eq.(jl]): 

^ = ^^1^E^X^'(^), (19) 

n 

using experimental values (T^^^(s). But comparing (fT9!) with definition (fTOj) . 

nln2~-ln<^'P(s), (20) 

it seems that last one gives more definite value of z^^P{n, s) than the "integral" equation 
(TT9l) especially since the statistical errors are large in the VHM region and the theoretical 
correction to Eq. (!20l) are small, ~ 1/n. 

Summarizing the results we conclude: if the energy is sufficiently large, i.e. if Zmax is 
sufficiently close to z^ = 1, if the multiplicity is sufficiently large, so that ( (77| j is satisfied 
and z{n, s) can be sufficiently close to Zc, then one may have confident answer on the 
question: observable or not the first order phase transition in hadron/ion collisions. The 
heavy ion collisions are favorable to observe the phase transition. 
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